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Abstract. By using ab initio methods on different levels we study the 
magnetic ground state of (finite) atomic wires deposited on metallic surfaces. 
A phenomenological model based on symmetry arguments suggests that the 
magnetization of a ferromagnetic wire is aligned either normal to the wire and, 
generally, tilted with respect to the surface normal or parallel to the wire. From 
a first principles point of view, this simple model can be best related to the so- 
called magnetic force theorem calculations being often used to explore magnetic 
anisotropy energies of bulk and surface systems. The second theoretical approach 
we use to search for the canted magnetic ground state is first principles adiabatic 
spin dynamics extended to the case of fully relativistic electron scattering. First, 
for the case of two adjacent Fe atoms an a Cu(lll) surface we demonstrate that 
the reduction of the surface symmetry can indeed lead to canted magnetism. The 
anisotropy constants and consequently the ground state magnetization direction 
are very sensitive to the position of the dimer with respect to the surface. We 
also performed calculations for a seven-atom Co chain placed along a step edge 
of a Pt(lll) surface. As far as the ground state spin orientation is concerned 
we obtain excellent agreement with experiment. Moreover, the magnetic ground 
state turns out to be slightly noncollinear. 



PACS numbers: 75.10.Lp, 75.30.Gw, 75.75.+a 
1. Introduction 

Magnetic devices on the atomic scale became recently a subject of intensive 
experimental and theoretical research (see, e.g., the "viewpoint" drawn by Kiibler £Q). 
Understanding and design of the relevant physical properties - magnetic moments, 
magnetic anisotropy energies, thermal stability, switching - of atomic scaled magnets 
demand a detailed knowledge of their electronic and magnetic structure. For this 
reason a considerable amount of theoretical work has been published to investigate 
the - mostly noncollinear - magnetic ground state of free and supported metallic 
clusters H El El 13 El 
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Quite recently, Gambardella et al. reported well characterized experiments on 
linear chains of about 80 Co atoms located at a step edge of a Pt(lll) surface 
terrace UJ. At 45 K the formation of ferromagnetic spin-blocks of about 15 atoms 
was found with an easy magnetization axis normal to the chain and pointing along a 
direction of 43° towards the step edge. Stimulated mainly by this experiment, in the 
present work we present a study of the magnetic ground state of linear atomic chains 
deposited on a fcc(lll) host surface. We first focus on the origin of the unusual canted 
magnetism. In the case of two adjacent Fe atoms placed into a Cu(lll) surface we 
investigate how the orientation of the magnetization depends on the position of the 
dimer with respect to the surface. Then we calculate and analyze in some detail the 
ground state spin configuration of a finite Co chain at a Pt(lll) step edge. 

2. Theoretical methods 

Let us write the magnetic moment of a ferromagnetic system in terms of spherical 
coordinates, M= M (sin. 9 cos (f), sin sin 0, cos#), where < 9 < tt and < cf) < 2tt 
are the azimuthal and polar angles in a usual rectangular reference of frame, and 
assume that the magnitude of the magnetic moment, M , is independent on the the 
orientation. In case of a linear chain of atoms deposited along the x axis of an fcc(lll) 
lattice (see, e.g., figure ^| the system has at best one symmetry operation, namely, a 
mirror symmetry with respect to the (y, z) plane, therefore, invariance of the energy 
implies up to second order in the magnetization that 

E(9,4>) = E Q + if 2 ,i cos20 + #2,2(1 - cos 29) cos 2$ + K 2 ,3 sin 29 sin (j) , (1) 

where #2.1 (i = 1,2,3) are so-called anisotropy constants. Solving the corresponding 
Euler-Lagrange equations results that, depending on the actual values of the 
anisotropy constants, the easy magnetization axis corresponds either to <p = tt /2 and 
9 E {9q,9 q + tt/2, 9 Q + tt}, where 

or to 9 — tt/2 and <fi — 0. Clearly, only in the special case of #2,3 = can the ground 
state magnetization point along the z axis (perpendicular to the planes). 

The so-called magnetic force theorem (MFT) represents a straightforward and 
relatively simple way to calculate anisotropy constants as based on the local spin 
density approximation (LSDA). Here, a self-consistent calculation is carried out for 
only one selected orientation of the magnetization. Then, by keeping these potentials 
and effective fields fixed, the orientation of the spin-magnetization (in LSDA parallel 
to the effective field) is varied, whereby - neglecting further self-consistency - only 
the single-site (band) energy is considered, see reference ^HI for several applications 
of this method in ordered and disordered layered systems. Note that, in principle, also 
the magnetic dipole-dipole interaction energy has to be added to the bandenergy. For 
the case of small nanostructures the estimated magnitude of this energy is, however, 
by at least one order less than that of the bandenergy. 

A numerically efficient tool to search for an equilibrium spin arrangement is to 
trace the time evolution of the spin moments until a stationary state is achieved. 
According to the so-called first principles adiabatic spin-dynamics (SD) founded by 
Antropov et al. 11 , for a system with well-defined local (atomic) moments the 
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evolution of the time dependent orientational configuration, {ej(t)}, is described by a 
microscopic, quasi-classical equation of motion, 



x(e,xBf) , (3) 



where B^^ is an effective magnetic field averaged over cell i, 7 is the gyromagnetic 
ratio and A is a damping (Gilbert) parameter. Following the arguments of Stocks 
et al. |12l ll.'ij in this equation at any moment of time the orientational state has 
to be evaluated within a constrained density functional theory (DFT). Here a local 
constraining field, B™" ensures the stability of a non-equilibrium orientational state. 
This implies that the internal effective field that rotates the spins in the absence of 
a constraint and, therefore, has to be used in equation © is just the opposite of the 
constraining field ^2- By merging with the locally selfconsistent multiple scattering 
(LSMS) method SD has been applied so far to bulk metals and alloys [121 1131 HI] and, 
very recently, to interfaces |15| . 

In order to deal with exchange splitting and relativistic scattering on equal 
theoretical footing we combined the first principles SD scheme based on constrained 
DFT by solving the Kohn-Sham-Dirac equation, 

[cap + f3mc 2 + V(r) + fi B P <r-(B xc {r) + B co "(r)) - E] = , (4) 

where a and /3 are the usual Dirac matrices, <x are the Pauli matrices, V(r) stands 
for the Hartree and the exchange-correlation potential, while within the local spin 
density approximation (LSD A) B a;c (r) is an exchange field interacting only with the 
spin of the electron. Equations @ and (0J) form the very basis of a relativistic spin- 
only dynamics, inasmuch no attempt is made to explicitly trace the time evolution of 
the orbital moments. 

In conjunction with both the MFT and the SD we applied the multiple scattering 
Green's function embedded cluster method developed by Lazarovits et al. 16 . In 
here, first a self-consistent calculation is carried out for the surface system in terms 
of the relativistic Screened Korringa-Kohn-Rostoker (SKKR) method |17|. and then 
the nanostructure is embedded into this host according to the equation 

7*(e) = 7-(e) [I - (t^e)" 1 - t^e)" 1 ) -f(e)] 1 , (5) 

where T r (t r ) and T°(t c ) are site-angular momentum matrices of the scattering 
path operators (single-site t operators) of the host surface system and the cluster, 
respectively, and e is the energy. By solving also the corresponding Poisson equation 
with appropriate boundary conditions a selfconsistent calculation for the selected 
cluster can be performed that takes full account of the environment ^Hl- It is 
important to underline that this description does not rely on periodic boundary 
conditions applied to the (embedded) cluster. In all calculations we used the atomic 
sphere approximation (ASA) and an angular momentum expansion up to £ m ax = 2. 



3. Two Fe impurities at Cu(lll) surfaces 

The simplest system that reduces the symmetry of the fcc(lll) layers to only one 
mirror plane is a pair of nearest neighbour defects as illustrated in figure ^ We 
embedded two such Fe impurities on top, i.e., into the first vacuum layer (labelled by 
S + 1), into the surface (S) layer and into the subsurface (S — 1) layer of a Cu(lll) 
surface. First a self-consistent SKKR calculation has been performed by relaxing the 
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potentials of 6 Cu layers and 3 empty sphere layers in order to describe the surface 
region from the bulk to the vacuum. Then, by employing equation JSJ, for the three 
above cases the Fe dimers were calculated self-consistently. Since these calculations 
are intended provide only with qualitative predictions, for simplicity, we neglected 
relaxations of the potentials of the host atoms, 
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Figure 1. Schematic top view of two impurities (big shaded circles) placed into 
fcc(lll) layers as nearest neighbours. The host atoms in the same layer are 
displayed by small shaded circles, those in the layer above by full circles, while 
those in the layer below by empty circles. 



As mentioned before, the present implementation of our SD scheme serves (only) 
for finding the magnetic ground state of the system, therefore, it is sufficient to consider 
only the second (damping) term on the right hand side of equation }3"|. The evolution 
of the spin orientation is then measured on a time scale with a unit (time step) of 1/A. 
In figure [21 the evolution of the 9 and 4> angles and the magnitude of the constraining 
field are plotted in this artificial time scale for one of the Fe atoms. Note that during 
the SD procedure, the magnetic configuration of the two Fe atoms was confined to 
be symmetric with respect to the (y, z) plane. (We checked, however, that the final 
result is independent on the starting configuration.) 
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Figure 2. Evolution of the direction of magnetization and the magnitude of the 
constraining field according to a spin-dynamics calculation for two Fe impurities 
placed into different layers of a Cu(lll) surface: ■ S + 1, A S, O 5 — 1. For 
better visibility, the corresponding results for the first 300—1000 iterations are not 
shown. Displayed are only the data for one of the Fe atoms (see text). 



Actually in a few time steps, for all the three layer positions the magnetic state 
of the two Fe atoms became nearly ferromagnetic and perpendicular to the line 
connecting the two impurities (<j> — 90°). A satisfactory convergence was, however, 
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achieved only after thousands of time steps later, when - as can be seen from figure[21- 
the constraining fields converged to zero. The final magnetic states can be summarized 
as follows: 6 = 0.73°, 66.8° and 13.1°, as well as <f> = 89.8°, 89.9° and 89.9° for one of 
the Fe atoms in layers 5+1, 5 and 5 — 1, respectively, and a symmetric orientation for 
the other Fe atom. Clearly, these magnetic ground states are in qualitative agreement 
with the predictions of the simple phenomenological theory, see equations and J2J . 
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Figure 3. Bandonergy differences from MFT calculations for two Fe impurities 
at a Cu(f If) surface: ■ S + f , A 5, O S - f ; a) 4> = 90° , b) <j> = 0. Solid lines 
show curves fitted to equation with the parameters contained by table HI 



Table 1. Anisotropy parameters (in units of meV), see equation Q, fitted to 
energies calculated within the MFT for two Fe impurities placed at different layers 
of a Cu(ff f) surface. The azimuthal angle corresponding to the global minimum 
of the energy is also displayed. 
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K 2 ,3 
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S+l 


-0.5f 


-0.39 


-0.0026 


0.083° 


S 


-0.f8 


0.23 


-0.079 


62.0° 


S-l 


-0.21 


O.fO 


-0.065 


f5.7° 



We carried out MFT calculations by using the output of the SD calculations. 
Figure [3] shows (band)energy curves when scanning the (uniform) direction of the 
magnetization through the paths < 9 < 180° for </> = 90°, i.e., in the (y, z) plane 
and for = 0°, i.e., in the (x, z) plane. All the calculated data fit almost precisely the 
function in equation Q with the parameters listed in tabled For all the three cases 
the minimum of the energy is found in the (y, z) (symmetry) plane at an azimuthal 
angle, 9, also shown in table ^ These angles coincide remarkably well with those 
obtained from the SD calculations. It should be noticed, however, that by using the 
output of self-consistent calculations with a magnetization fixed along the z axis, in 
particular, for an Fe dimer placed into layer 5 we obtained an apparently different 
ground state orientation (9 — 46.5°). This clearly demonstrates that the applicability 
of MFT is quite limited. 

It is obvious that a pair of Fe atoms in layer 5 + 1 represents the case of strong 
perpendicular anisotropy with a corresponding anisotropy energy of about 0.9 meV/Fc 
atom. This clearly explains the relatively fast convergence of the SD scheme, while 
in the other two cases due to the smaller spin-orbit coupling the convergence was 
much slower (see figure |2J). For Fe dimers immersed into layers 5 and 5 — 1 a change 
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of the orientation from the z to the x direction costs much larger energy than a 
corresponding variation within the symmetry plane, because the otherwise remarkably 
reduced anisotropy constants, K^i and K 2 ,2, differ in sign. Concomitantly, the relative 
increase of the magnitude of if 2, 3 gives rise to a canted ground state, see equation J2J). 
It should be noted that, similar to as discussed in Ref. [T^], the change of the anisotropy 
constants, K 2 ,i with respect to the position of the Fe dimer can be related to the 
different hybridization between the electronic states of the Fe and the Cu atoms. 

4. Finite Co wire at the step edge of a Pt(lll) surface 

We first performed a calculation treating 8 layers of Pt selfconsistently together with 4 
layers of vacuum. Then, a seven-atom chain of Co together with 10 empty (vacuum) 
spheres were embedded into the topmost Pt layer as schematically indicated in figure 
0] Simultaneously all the nearest neighbours of the Co atoms were re-embedded into 
the respective Pt or vacuum layers to allow for the relaxation of potentials around 
the Co chain. Therefore, an embedded cluster of a total of 55 atoms was treated 
selfconsistently. 



Figure 4. Schematic view of the geometry of a seven— atom Co chain along a 
Pt(lll) step edge. Full circles: Co atoms, shaded circles: Pt atoms, open circles: 
empty spheres, a): side view, b): top view of the surface Pt layer with the Co 
chain. The cluster embedded is indicated by solid lines. The coordinate system 
gives reference to the azimuthal and polar angles, 9 and 0, that characterize the 
orientation of the magnetization. (Note that in Ref. |2] a different coordinate 
system and the opposite notation for the angles is used.) 



For each Co atom in the chain, in figure [5] the evolution of the 8 and <j> angles is 
plotted for the first 100 time steps in the artificial time scale mentioned in section [21 
Initially the directions of the atomic magnetic moments were set by a random number 
generator. It can be seen that after some oscillations both the 9 and the <f> angles for 
all the Co atoms quickly approach to nearly the same value for all the Co atoms, i.e., 
as in the previous case of Fe dimers at a Cu(lll) surface, to a nearly ferromagnetic 
configuration. The initial rapid oscillations that can be observed in figure El are the 
consequence of the relatively large constraining fields caused by large exchange energies 
whenever the moments point into very different directions. In about 1000 time steps 
the (j) angles, converged to 90°, with a spread of less then 1°, i.e., normal to the chain 
and the 9 angles converged to nearly 42°. All these results are in excellent agreement 
with experiment [§]. 

We also performed MFT calculations with the selfconsistent potentials and fields 
obtained from the SD procedure. The fitted anisotropy parameters K 2t i = —0.16 
meV, i^2,2 = —1.06 meV, and -^2,3 = —4.81 meV not only fairly well reproduce the 
easy axis, 8 = 38° and = 0°, but also result in a value of the anisotropy energy 
(defined as the energy difference between the hard and the easy axes) of 1.42 meV/Co 
atom that satisfactorily compares with that derived from experiment (2.0 meV) [Jj]. 
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Figure 5. Evolution of the angles 8, part a), and <j>, part b), defining the 
orientation of the spin moments for the seven Co atoms in the finite chain depicted 
in figure [1] The symbols refer to the following Co atoms numbered from the left 
to the right in part b) of figure^] ■ 1, O 2, ▲ 3, ♦ 4, A 5, • 6, □ 7. Shown 
are only the first 100 time steps. 



Table 2. Calculated magnitudes and orientations of the spin and orbital moments 
in a seven-atom Co chain along a Pt(lll) step edge. 



atom Spin moment Orbital moment 

moment (/zb) 0(deg) moment (/is) ©(deg) 



1 


2.23 


41.1 


0.25 


39.1 


2 


2.18 


42.5 


0.20 


41.5 


3 


2.18 


42.3 


0.19 


40.1 


4 


2.18 


42.4 


0.20 


41.3 


5 


2.18 


42.3 


0.19 


40.2 


6 


2.18 


42.5 


0.20 


41.5 


7 


2.23 


41.1 


0.25 


39.1 



Extracted from the final (equilibrium) state, the size and the azimuthal angle 9 
of the spin and orbital moments for each Co atom are shown in tabic While the 
calculated spin moments for the inner Co atoms (2.18 hb) are in good agreement 
with the value deduced from experiment (2.12 /zb) and also with other theoretical 
studies on infinite wires [^||2Jj], the edge of the wire is characterized by larger spin (and 
orbital) moments [2U- Although our calculated orbital moments for the inner atoms 
(0.19-0.20 /j,b) are larger than the corresponding values from other LSDA calculations 
(0.16 fis ,E| and 0.15 hb |20|)> they are still much too small when compared to the 
experimental value (0.68 /is) 0. Note that including orbital polarization scheme or 
using the LDA+U method a value of 0.92 [is and 0.45 /ig [201 can be obtained. 

As can be inferred from table the spin moments of the inner atoms are fairly 
parallel, those at the end of the chain however, are off by more than 1°. This can be 
associated with the anisotropy energy contributions being larger at end of the chains 
than inside as found for finite Co wires deposited on a Pt(lll) surface [22- It can 
also be seen in table |21 that the orbital moments oscillate stronger in magnitude and 
orientation than the spin moments. As pointed out by Jansen |22| . within the DFT 
the spin and orbital moments are required to align only when the ground state refers 
to a high-symmetry direction. This is, however, not the case for the Co wire since the 
ground state orientation is not directly induced by symmetry. 
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5. Summary 

In this work we presented calculations of the magnetic ground state of linear atomic 
chains placed onto surfaces in terms of the magnetic force theorem and an ab-inito 
spin dynamics scheme. We found that due to the variation of the anisotropy constants 
the canted magnetic state obtained for a Fe dimer at a Cu(lll) sensitively depends 
on its position. In excellent quantitative agreement with experiment, we obtained a 
canted ground state for a finite Co wire along a Pt(lll) surface step edge. We also 
found that this magnetic state is noncollinear: a feature that is expected to play a key 
role in nanostructures of more complex geometry. 
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